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Abstract 

The topic of this review is the current state of our knowledge about the degree 
of initial equilibration between electrons, protons and ions at supernova remnant 
(SNR) shocks. Specifically, the question has been raised as to whether there is an 
inverse relationship between the shock velocity and the equilibration similar to the 
relationship between equilibration and Alfven Mach number seen in interplanetary 
shocks (Schwartz et ah, 1988). This review aims to compile every method that has 
been used to measure the equilibration and every SNR on which they have been 
tested. I review each method, its problems and uncertainties and how those would 
effect the degree of equilibration (or velocity) inferred. The final compilation of 
observed electron to proton temperature ratios as a function of shock velocity gives 
an accurate, conservative picture of the state of our knowledge and the avenues we 
need to pursue to make progress in our understanding of the relation between the 
velocity of a shock and the degree of equilibration. 
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1 Introduction 


The nature of electron and ion heating behind collisionless shocks remains an 
open question in shock physics. While the overall kinetics of a shock can be 
described simply in terms of the Rankine- Hugoniot solutions to the equa tions 


for conservation and continuity (see e.g. McKee and Hollenbach . IQSflll . the 
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fraction of the shock’s kinetic energy that is transferred to the thermal and 
cosmic-ray populations of electrons and ions is unknown. Theoretically, any¬ 
thing between (Te/Tp)o ~ rrie/mp to rapid full equilibration has been proposed. 
However mounting observational evidence points to an intermediate degree of 
initial electron heating, where the disparity between electron and ion heating 
increases with shock velocity. This review examines these measurements in 
order to point the way towards further progress. 


Modeling non-linear acceleration of particles in high Mach number collisionless 
shocks remains a challenge. There are numerous analytically derived instabil¬ 
ities that could be important, but it is not necessarily straight-forward to 
interpret particle simulations in terms of these processes. While comparing 
the simulations to observations of solar wind shocks has been quite successful, 
the results may not have a strong bearing on the important processes in super¬ 
nova remnant (SNR) shocks which are signihcantly faster and less dense. For 
a r ecent, thorough review of the open problems in collisionless-shock physics 


see 


Lembege et ahl (|2nn4 h For the purpose of this review, it is important to 


keep in mind that the electro-magnetic waves that dehne these instabilities 
and the shock itself will affect the ions and electrons differently due to the 
large inertial difference. 


Given the current ambiguity regarding the partitioning of energy at a shock, 
interpretations of SNR observations commonly assume purely thermal particle 
distributions and consider the extreme cases of minimal equilibration (mass- 
proportional heating) or full equilibration (all species rapidly heated to the 
same temperature). ^ If cosmic ray production is included, the case consid¬ 
ered is usually one with an extremely high injection efficiency of a thermal 
seed population into the cosmic ray acceleration process so that the cosmic 
ray flux at Earth can be explained. For observations that only measure one or 
two populations, the choice of equilibration can drastically affect the inferred 
age, explosion energy or distance. Furthermore, if equilibration is a function 
of shock speed, then even spanning both scenarios may be misleading because 
the evolution of the initial shock equilibration will effect the dynamics of the 
remnant (for instance consider the internal pressure). Calibrating a relation¬ 
ship between shock velocity and the efficiency of heating different particles to 
thermal and cosmic ray energies would be a fantastic boon for SNR studies if 
a clean relationship exists. 


There are three generic methods for determining the temperature of parti¬ 
cles behind a shock: (1) bremsstrahlung continuum emission, (2) the thermal 
broadening of a line from a particular species, and (3) the flux ratios of lines of 
either a single element or between elements. Each of these methods has its own 


Bocchino et al. mm) turned this argument around to estimate a degree of equi¬ 


libration based on given Sedov parameters for the Vela supernova remnant 


2 









pitfalls. (1) Both thermal bremsstrahlung and synchrotron radiation (from a 
high energy population of electrons) produce continuum emission in the X- 
rays, and their relative contributions are not always possible to disentangle. 
Furthermore the bremsstrahlung continuum in the X-rays samples emission 
from an electron population signihcantly downstream from the shock. Assum¬ 
ing the electrons and protons were out of equilibrium the temperatures will 
equilibrate downstream on the timescale of coulomb collisions. The initial elec¬ 
tron temperature must then be inferred using the ionization state to calibrate 
the coulomb timescale post-shock. (2) The line broadening is a function not 
only of the temperature of that species but also the velocity structure of the 
shock if the shocked hlament is not resolved. (3) The temperature dependence 
of ionization and excitation by electrons and protons can be used to diagnose 
Tg or Tp from the line flux ratios. The fluxes also depend on the ion fractions 
or relative abundances which must be assumed if they cannot be derived. We 
will examine many manifestations of these generic methods and problems as 
the observational evidence unfolds. 

The body of this review is organized by supernova in order of increasing ve¬ 
locity to emphasize the necessity of measurements in multiple wave-bands to 
corroborate the degree of equilibration. Although different authors use dif¬ 
ferent parameters to express the degree of equilibration, the ratio of the ini¬ 
tial electron to proton temperatures at the shock front (Tg/Tp) will be used 
throughout this review. These studies all concentrated on the equilibration 
at the outer blast wave of each remnant. Further complications must be con¬ 
sidered with regard to the energetics of a reverse shock into the pure metal 
ejecta, where the energetics are not dominated by the protons, cold electrons 
are continually released by ionization, and turbulence may generate a high 
magnetic held. 


2 The Cygnus Loop and RCW 86 


T he ratio of broad-to - narro w component hux in the Ha and H/? lines was used 
bv iGhavamian et ahl ()2nniri to measure the electron-proton temperature ratio 
behind shocks in the Cygnus Loop and RCW 86. The t heory behind th e broa. d 
and narrow Balmer line emission was hrst laid out by Chevalier et ah ( 198flll . 
Directly behind a collisionless shock there can exist a population of cold neu¬ 
tral ions that were not affected by the shock passage. These neutrals will emit 
Ha and H/? Balmer lines if they are collisionally excited before being ion¬ 
ized. This produces a narrow component of emission, whose width rehects the 
temperature of the pre-shock gas. The cold neutrals can also charge exchange 
with the shock-heated protons resulting in a population of hot neutrals. A 
broad component to the Balmer lines arises from this hot neutral population, 
where the width of the broad line represents the temperature of the post-shock 
protons. 
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Fig. 1. The Ha and H/3 profiles of ( a) the Cygnus Lo o 
compact knot in Tycho’s SNR from lOhavamian et aT 


) p. (b) RCW 86 and (c) a 
(12001 i. Note that the Ha 


profile has higher signal to noise to constrain the broad-to-narrow ratio than does 
H/3, with the possible exception of the Cygnus Loop. 


Figures la and lb from iGhavarnian et abl (j2nnill present the Ha and H/3 pro¬ 
files from the Cygnus Loop and RCW 86. ^ The post-shock proton tempera¬ 
tures derived from the widths of the Ha broad components are Tp = (1.5 ± 
0.26) X 10® K for Cygnus and Tp = (6.9±0.2) x 10® K for RCW 86. To estimate 


the de gree of equilibration between electrons and protons iChavamian et al. 
(1200 ill modeled the broad and narrow emission including the temperature 
dependence of charge exchange and excitation of the fast and slow neutrals 
by electrons and protons. The predicted broad-to-narrow ratios are plotted 
for different values of t he initial neutral fraction in hgures 2a and 2b from 
Chavamian et a].l (j200111 . Using the H/3 ratio for Cygnus and the Ha ratio for 
RCW 86, they report constraints of Te/Tp= 0.67 — 1.0 and Te/Tp= 0.25 — 0.43, 
respectively. Assuming cosmic ray production was unimportant to the energet¬ 
ics of the shock, these equilibrations and the broad component widths imply 
shock speeds of 300-400 km s~^ (Cygnus) and 580-660 km s~^ (RCW 86). 

However, it i s worth noting frq r n hgu re 2a that the Ha broad-to-narrow ratios 
predicted by Chavamian et al. ( 200 ill for the Cygnus loop lie 10% higher than 
those observed. This level of discrepancy is readily explained by contamination 
of the narrow component by emission from a shock precursor. Some fraction of 
the hot neutrals will escape ahead of the shock before being ionized. Further¬ 
more, in a cosmic-ray modihed shock, some fraction of the cosmic rays will also 
escape into the pre-shock gas. The hot neutrals and cosmic ray particles will 
ionize and excite the medium ahead of the shock. This non-linear effect was not 


included in the models of IChavamian et ahl ()2nnill and in fact, modeling the 
balance between the production of hot neutrals and the continued existence of 
neutral atoms at the shock front is a work in progress (K. Korreck, personal 
communication). However, it is clear that precursor emission will increase the 
narrow component flux, shifting all the models to lower broad-to-narrow ra- 


^ The data are presented without the model so that the reader can better evaluate 
the reliability of the broad component measurement for each remnant in both Ha 
and H/3. 
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Fig. 2. The predicte d Ha and H/3 broa d -to-na rrow ratios as a function of “fractional 
equilibration” from lOhavamian et all (12001 ) for velocities appropriate for (a) the 
Cygnus Loop, (b) RCW 86 and (c) Tycho. Note that the model predictions are 
inconsistent with the higher signal-to-noise Ha broad-to-narrow ratio for the Cygnus 
Loop and Tycho’s SNR. 


tios. Given this uncertainty it is unclear whether a measured broad-to-narrow 
ratio lower than all of the models provides any constraint on T^/Tp at all. In 
fact in sections 3 and 4 we will see examples of shocks in which the broad-to- 
narrow ratio is signihcantly lower than any of Ghavamian’s models, strongly 
indicating precursor emission. At present there is no reason to assume that 
this level of precursor contamination is missing in the Gygnus loop or RGW 86 
just because their broad-to-narrow ratios happen to marginally intersect the 
predicted curves. 


On the other hand, there is corroborating evidence for a high degree of equili¬ 
bration in the Gygnus Loop. Raymond et ah ( 200311 examined the far ultrav i- 
olet spectrum at the same position as was studied in Ghavamian et al. ( 200311 . 
The width of the O IV A1032 line gives an upper limit on the oxygen temper¬ 
ature of 2.7 X 10® K (any velocity structure broadening would decrease this 
es timate). This is wi t hin a factor of 2.5 of the hydrogen temperature measured 
Ghavamian et al. IjlEil). Given that mass-proportional heating would lead 


m 

to a factor of 16 difference between the temperatures of oxygen and hydrogen, 
this measurement indicates a high degree of equilibration between ions at the 
shock. This lends credibility to the claim of high electron-ion equilibration. 
However, while the electron heating mechanisms may be tied to the tempera¬ 
ture of the ions, there is still a three orders of magnitude difference in mass. 
Hence a high degree of equilibration between the ions themselves does not nec¬ 
essarily imply a similarly high degree of equilibration between the electrons 
and ions. 
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Fig. 3. Initial electron to proton temperature ratio as a function of s hock speed for 
five apertures in DEM L71. Published results ( Rakowski et al. . 20031 1 are plotted as 


solid crosses, new result as dashed. The new analysis incorporates lOOks more data 
which enables constraints on the abundances for each region. The offset between the 
old and new results illustrates how an assumption about the underlying abundance 
set can bias the measured temperature. 


3 SNR DEM L71 


A combination of optical and X-ray observations was required to constrain 
the electron-proton temperature ratio around th e blast-wave of SNR DEM L71 
(jOhavamian et ab , 20031: Rakowski et al. , 2003[l . This Large Magellanic Cloud 
(LMC) remnant has many non-radiative Ha hl aments around its outer rim, 


ranging in velocity fro m 500 — 1000 km s IChavamian et al.1 (|2003ll and 


iR.akowski et al. (2003) chose apertures along the blast-wave to combine areas 
of similar optical properties into large enough regions for meaningful spectra 
in a 45.4 ks Chandra observation. The optical measurements of the Ha broad 
component width (to obtain Tn) and attempts to model the broad-to-narrow 
ratio are presented in IChavamian et ( 2003ll . Without the addition of pre¬ 
cursor emission the predicted broad-to-narrow ratios never drop below 1.0, 
but the measured flux ratios are never more than 0.74 (generally around 0.5). 
Until full precursor modeling is available no conclusion can be made about 
the degree of equilibration purely from the optical data. 


In Rakowski et al. ( 20031) . we extracted Chandra X-ray spectra from three 


nested regions behind the blast-wave at each aperture to measure the evolu¬ 
tion of the electron temperature. Our hnal analysis of the initial ratio, Te/^p, 
utilized a pl anar shock model with t he initial proton temperature hxed to the 
values from Chavamian et al. ( 20031) . The electron temperature as a function 
of distance behind the shock could then be parameterized in terms of the initial 
degree of equilibration. The published results are shown as solid crosses in £g- 
ure 3. Except for the slowest shock (region 5), all of the shocks were consistent 
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with no equilibration in the preliminary tests. Combining X-ray and optical 
results allowed us to identify one shock (region 4) with Te/Tp~ 1/2 whose 
spectrum was inconsistent with either full or minimal equilibration given the 
measured proton temperature. Limited statistics dictated that the abundances 
and column density be held constant in the X-ray hts. The chosen average 
abundance set will introduce a bias given the inverse correlation between the 
abundance and temperature when both are allowed to vary. The errors shown 
include an estimate of this systematic error as well as the measurement un¬ 
certainty in Tp. An additional 100 ks of Chandra observing time has allowed 
us to eliminate this systematic uncertainty by measuring the abundances and 
column density at each aperture. These data conhrm that our shock model is 
sound, i.e. the temperatures and ionization timescales of the nested regions 
do indeed evolve when £t separately. The hrst of the updated measurements 
of Tg/Tp is shown by the dashed cross in figure 3 conhrming an intermediate 
degree of equilibration for region 4. While the uncertainty remains similar, the 
result itself is no longer biased by an arbitrary choice of abundances. 


4 Tycho’s SNR 


H a and H/3 emission from a compact knot in Tycho’s SNR were also included 


m 


Ghavamian et ah ( 200111 (see section 2). The spectra and model predic- 


tions for the broa d To-na rrow ratios are shown here in hgures Ic and 2c from 
Ghavamian et al.l (j200jJ). The Ha prohle has considerably better signal to 
noise than H/3, but the Ha broad-to-narrow ratio falls belo w predictions even 


after accounting for some diffuse emission above the knot. iGhavamian et ah 


( 200111 therefore use the less well-constrained H/3 ratio to estimate T^/Tp < 


0.1. The predicted broad-to-narrow ratios in hgure 2c certainly favor a low 
equilibration. The narrow component contribution from the cosmic-ray pre¬ 
cursor would have to be more than one third of the total narrow Ha emission 
in order to be consistent with a T^/Tp ratio greater than 0.2. However, a more 
precise statement on T^/Tp awaits a full precursor model to understand the 
difference between the predicted and observed Ha ratio. 


Hwang et al. ( 200211 studied the Tycho SNR blast-wave with Chandra. Spectra 
extracted from narrow regions along the outer rim exhibited little or no line 
emission. Assuming this emission is thermal in origin a very young ionization 
age is needed to explain the lack of line emission given a htted electron temper¬ 
ature around 2 keV. Given the low timescale, the electron temperature they 
measure should not have evolved since the initial shock heating. This electron 
temperature can be compared with a shock velocity derived (at a distance of 
2.3 kpc) from radial expansion measurements in the radio (3(300-40 00 km s“^ 


Revnoso et al. . 1997 1 and X-rays (4600 ± 400 km s ^ Hnghe4 12 OOOII to derive 


the ratio of the electron to mean temperatures. There remains a debate over 
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which expansion measnrement is representative of the present shock speed, 
however the electron-to-mean temperatnre ratio would be 0.10—0.20 assum¬ 
ing the mean radio expansion or 0.054—0.10 for the X-ray expansion speed. 
Combining both these ranges yields a limit on Tp/T^, of 0.03 —0.12. This is 


in good agreement with the results of iGhavamian et ahl (1200111 . however at a 
much faster shock within the same remnant. Hwang et ah I (120021) note that the 
low measured electron temperature does not require a significant cosmic ray 
population to explain the gap between the mean temperature for the measured 
shock speed and the electron temperature. However, neither does it exclude 
the possibility that a significant fraction of the shock energy has gone into 
accelerating cosmic rays. 

The above conclusions were derived under the assumption that the blast-wave 
emission is thermal. However, synchrotron emission arising from high energy 
electrons accelerated to cosmic-ray energies at the shock front is an equally 
natural explanation for such a featureless spectrum. The presence of hard X- 
ray emission up to 30 keV requires that a substai itial portion of the 0.5-10 
keV continuum be no nthermal. fsee sectio n 4.2 of Hwang et ah ( 20021) for a 
full discussion) In fact, Hwang et ah ( 20021) were able to successfully model the 
line-free rim spectra with a synchrotron emission model that required consis¬ 
tency with radio measurements. They found that if the emission is primarily 
non-thermal, any additional thermal component would have to have a con¬ 
siderably lower electron temperature than was found assuming a completely 
thermal origin. This would be consistent with a picture where a large portion 
of the shock energy has gone into accelerating cosmic rays, leaving behind a 
low temperature population of thermal electrons, which may or may not be 
equilibrated with the thermal population of ions. 


5 SN 1006, fully unequilibrated shocks 


SN 1006 is the one remnant where the equilibration between particles has been 
truly well-studied. Four lines of evidence all indicate a low degree of equili¬ 
bration or nearly mass-proportional heating. Any one estimate of Tl/Tp ma 


not be convincing, either because of its indirectness (jRavmond et ahl Il99 


I 


Vink et al.L 12003) or its dependence on arguments regarding H ie abundances 

However, the 


Ghavamian et ah . I 2 OO 2 I: Ihaming et al.L ^ 


or ion fractions 

convergence of four largely independent measurements is compelling. 


Raymond et ah ( 19951) observed the ultraviolet spectrum from a 197" long 
region along the north-west shock measuring emission lines from H I, He ll, C 
IV, N V, and O vi. They found that lines from all species show similar widths 
consistent with a speed of 230 0 km s ~^ indicating that ion-ion equilibration 
is ineffective. In Teaming et al. ( 19961) . they modeled the line ratios for this 
































































ultraviolet spectrum. From the ratio of two lines, one excited primarily by 
electrons, He ll A1640, and one excited by protons, C iv A1550, one can infer 
the ratio of electron to ion temperatures. Assuming solar abundances and a 
factor of two carbon depletion onto grains, this flux ratio indicates Te/Tion< 
0.05. However, since the flux ratio is obviously dependent on the assumed 
relative abundances, the authors admit that a ratio of Te/Tion= 0.2 cannot be 
ruled out. 


Ghavamian et ah ( 2002[l re-examined the optical spectrum of the northwestern 


shock. The most useful diagnostic pair of lines were He l A6678 and He ll A4686. 
The ratio of these two lines depends on the initial neutral fraction and the 
ratio of the He l excitation to ionization rates. The excitation to ionization 
rate drops sharply above Te > 10® K, implying that a low electron temperature 
is needed to produce measurable He l emission. A high neutral fraction is also 
necessary to replenish the supply of He l which is still ionizing faster than 
it emits photons. Quantitatively, comparing the electron temperature derived 
from modeling the He l to He ll ratio with the proton temperature measured 
with the Ha broad component width, yields a limit of Te/Tp< 0.03. The 
other measured line ratios, H l/Ha and the broad-to-narrow ratios of Ha and 
H/d, while less directly diagnostic, are all consistent with a low equilibration, 
Te/Tp< 0.07. 


I Vink et ( 2003ll studied the X-ray grating spectrum of a compact knot in 
the northwest using XMM-Newton. The broadening of a complex of O v, O 
VI and O VII lines imply an oxygen temperature of kT. 


. . .... . oxyp-en _ _ 

If there has been no ion equilibration, as suggested bv iRavmond et ahl ( 199511 . 
this oxygen temperature translates into a shock velocity of Vg > 4000 km 
s“^. If ion equili bration has occurr ed that would require an even faster, more 


= 530 ± 150 keV. 


energetic shock. Vink et ah ( 2003ll use a CCD-resolution X-ray spectrum, ex¬ 
tracted from the same XMM-Newton observation, to measure an electron 
temperature of kT^ ~ 1.5 keV. The ratio of the electron and oxygen temper¬ 
atures imply an electron-to-proton temperature ratio of ~0.05 for the case 
of no ion equilibration, or even smaller values for a shock where the ions are 
equilibrated. 


All of the measurements to date of SN 1006 indicate that the electrons and 
ions are far from being in temperature equilibrium. However, one of the most 
remarkable things about SN 1006 is the incomplete shell of bright radio and 
X-ray synchrotron emission in the northeast and southwest indicating sites of 
efficient cosmic ray acceleration. All measurements concerning equilibration 
have been made for the northwestern rim, which does not seem to emit X-ray 
synchrotron emission. The situation for the synchrotron rims may be quite 
different. 
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Fig. 4. Combin ed X-ray line profile from the grating spectrum of SN 1987A, re¬ 


produced from iMichael et al.l (|2nn2f) figure 8. Single (solid) and double (dashed) 


Gaussian fits are over-plotted. 


6 SN 1987A 


A single Chandra giaiing spectrum from an unresolved SNR, such as SN 1987A 
has the potential to measure the ion-ion and electron-ion equilibration simul¬ 
taneously. The line ratios of individual elements can be used to constrain the 
electron temperature, while the line widths are a function of the temperature of 
that element and the prohle of its bulk motion. The Chandra hig h energy trans¬ 


missio n grating (HETG) spectrum of SN 1987A, presented by iMichael et ah 


(|2002ll . displays Ka and La lines of N, O, Ne, Mg, and Si. Constraints on 
kTe and t he ionization t i mesca le from line ratios of individual elements were 
shown in Michael et ah ( 200211 . Unfortunately the grating spectrum is not 
deep enough to constrain the temperature and timescale fully, however the 
line ratios agree with CCD spectrum best-£t values of kT^ ~ 2.6 keV and 
rigt ~ 6 X 10^° cm“^ s. 

Ideally one would like to constrain each of the line widths separately, but 
again the statistics do not allow this. However, given that results from SN 
1006 indicate a low degree of ion equilibration at such fast shocks, it is not 
unreasonable t o assu me that all the ions have the same velocity. Therefore, 
Michael et ah ( 2002[l co-added the velocity prohles of the brightest lines of N, 
O, Ne, Mg, and Si to measure an “average” line-width. The combined X-ray 
line profile is reprinted here in hgure 4. Since SN1987A is not resolved by the 
gratin g, this line prohle inc ludes both thermal velocity and bulk motion. Fig¬ 


ure 


5 ( Michael et ah . 2002ll illustrates the line prohles produced by diherent 


geometries of the X-ray emitting plasma. Full equilibration between the ions 
produces sharper line prohles than no equilibration because of the greater bulk 
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Fig. 5. Line profiles predicted for differe nt geometries and de grees of equilibration 
for SN 1987A, reproduced from figure 9 of Michael et al. (20^). Solid lines show the 
bulk motion profiles. The dotted lines include the thermal broadening assuming full 
ion-ion equilibration. The dashed lines assume minimal equilibration at the shock 
front. 

motion contribution . The smooth, no eq uilibration cases are more similar to 
the line profile seen. Michael et al.l ( 200211 note that equatorial geometries are 
favored by other observations. For no ion equilibration and equatorial geome¬ 
tries Vg = 3400 ± 700 km s“^, consistent with radio expansion measurements. 


Similarly to what was done in DEM L71 by Rakowski et al. ( 2003ll . the velocity 
of the shock in SN 1987A can be used to constr ain the CCD spectrum in terms 
of the initial electron to ion temperature ratio. Michael et al. 1 (l2002ll hnd that 
if Vs ~ 3500 km s“^, then the CCD spectrum is ht best by T^/Ts = 0.11 ±0.02 
(where Ts is the averag e shoc k temperature, equivalently Te/Tp= 0.07±0.01). 
Unlike iRakowski et al. ( 2003ll . these constraints assume a hxed velocity and do 
not account for the range in velocities allowed by their line width calculations. 
However, the authors do note that hts with a very low electron temperature 
are inconsistent with the data, thus excluding the case of zero election-ion 
equilibration. 


7 IE 0102-72 


Up until this section, all discussions have assumed that the acceleration of 
cosmic rays is not an important factor in the shock dynamics. However it is 
widely believed that SNR shocks are the main source of the bulk of the cosmic 
ray distribution. If cosmic ray acceleration is efficient a signihcant portion of 
the shock energy may be diverted from the thermal populations. This should 
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have important consequences on the shock itself ( Decourcelle et ah . 200011 . 


The shock compression will rise due to the lower adiabatic index of relativis¬ 
tic particles, the forward and reverse shocks will be closer together and the 
temperature and density gradients will be more severe. As cosmic rays escape 
the shock they carry away shock energy, making these dynamic effects more 
dramatic. Neglecting the energetics of cosmic ray acceleration when consider¬ 
ing the equililbration of the thermal populations can give misleading results. 
If cosmic ray acceleration is important then results that are based on an ex¬ 
pansion velocity and an electron temperature may in fact have a much higher 
degree of equilibration bet ween the thermal populations (i.e. the Tycho re¬ 
sults of Hwang et al. . 200211 . Furthermore, inferring a velocity purely from the 
electron and ion temperatures will underestimate the shock speed if much of 
the shock energy is in relativistic particles. 


Hughes et ahl (|2000ll used Chandra calibration observations of IE 0102—72 to 
measure the expansion of the SNR since the earlier Einstein and ROSAT im¬ 
ages. With the known 50 kpc distance to the LMC, the 0.1% expansion per 
year implies a ~6000 km s“^ shock speed for the outer blast-wave. An X-ray 
spectrum of the blast-wave, extracted from the same observation, exhibits an 
electron temperature of less than 1 keV, 25 times smaller than the minimum 
average temperature implied by the measured shock velocity. There is no ini¬ 
tial Tg/Tp ratio for which the gulf between this average temperature and the 
electron temperature can be explained. Even after accounting for adiabatic 
decompression and only coulomb heating post-shock (up to the ionization 
timescale seen in the spectrum) the minimum electron temperature for such a 
fast shock is no lower than ~2.5 keV for the extreme case of mass-proportional 
heating. The discrepancy between the low electron temperature and high en¬ 
ergy of the shock requires the existence of another population into which the 
shock energy is flowing. If a signihcant fraction of the shock energy has gone 
into accelerating cosmic rays then our comparison of the electron tempera¬ 
ture with the shock velocity tells us nothin g about Te/Tp. In fact, non-linear 
models of shock acceleration ( Ellison . 200011 . at Mach numbers appropriate for 
IE 0102—72, exhibit electron temperatures near 1 keV and fully equilibrated 
electron and proton thermal populations. 


8 Conclusions 


The preceding results are compiled in hgure 6. Note the new constraint for re¬ 
gion 4 in DEM L71, obtained using a Chandra observation three times as long 
as the original. This region still requires an intermediate degree of equilibration 
even without assuming a particular abundance set. Recall a lso that corrobo¬ 
ratin g ion-ion equilibration measurements exist for SN in06(lRavmond et ah . 
I 995 L mass-proportional heating) and the Cygnus Loop ( Raymond et ah . 2003, 


12 







































0 1000 2000 3000 4000 5000 

shock velocity (km/s) 


Fig. 6. Compilation of current electron-ion equilibration measurements at SNR 
shocks. The solid boxe s for the Cygnus Loop, RCW 86 and Tycho, are from 
Ghavamian et al. ( 20011 1 modeling the broad-to-narrow ratios of Ha and H/3. The 


dashed lines are my estimate using only the H/3 measurement, however this esti¬ 
mate does not include the potentially significant effect of precursor contamination 
to the narrow c ompo nent. The grey box marked “Vela” shows the estimate by 


Bocchino et al.l (119901 1 of the current Tg/Tp, which did not take into account the 


amount of coulo mb equilibration that has already occurred. Crosses show DEM 
L71 results from Rakowski et al. (20^) using the broad Ha component to fix the 
proton temperature in an X-ray fit of Tg/Tp. The grey cross indicates the updated 
result from deeper observations which allow us to simultaneously constrain the el¬ 
emental abundances rather than fixing them to some average abundance set. The 
two SN 1006 resul ts indicate the most con straining of the optical or UV line ra¬ 
tio measurements ( Ghavamian et all . 2003[ l. and th e X-ray grating m easurement 
of the oxygen line width and electron temperature dViuk et all 200^1. T he solid 
box for SN 1987A shows the reported results from iMichael et al.l (l2002h using a 
fixed shock velocity to constrain Tg/Tpfrom the X-ray CCD spectrum, the dashed 
box gives an estimate which incl udes the error in their measurement of the shock 
velocity. The Hwang et al. I (l2002ll measurement plotted here, converted from their 
electron-to-mean temperature ratio for Tycho’s outer rim, compares the electron 
temperature measurement with both the radio and X-ray expansion velocities. For 
IE 0102—72, the low electron temperature could not be explained by a minimal 
degree of electron-ion equilibration at the shock front but rather required that a 
significant portion of the shock energy be going into the acceleration of cosmic rays. 
If the injection efficiency into the acceleration process is high, models predict that 
the remaining thermal populations may be equilibrated. 
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f 19991) of Tg/Tp between 0.4 and 1 for an approximately 600 km shock in 
the Vela SNR using ROSAT PSPC spectra. This was estimated by comparing 
the electron temperature of the hotter of two collisional ionization equilib¬ 
rium components with a proton temperature derived from their best estimates 
for the Sedov parameters for Vela. However a measurement of the ionization 
timescale is needed to calibrate the extent to which coulomb collisions have 
already equilibrated the electron and proton temperatures as was done for 
DEM L71 and IE 0102-72. 


< 2.5 X Tp). Also included in this hgure is the report bv iBocchino et ah 


An inverse correlation between shock velocity and efficiency of electron heat¬ 
ing has been suggested based on different subsets of these results. However, 
even if we accept that the high equilibration result for the Cygnus Loop is 
not an artifact of neglecting the pre-cursor emission, we are still left with one 
secure identihcation of an intermediate equilibration in the interesting range 
from 400 to 2000 km s“^ where the transition from near full equilibration to 
low equilibration should occur. Furthermore, once cosmic ray acceleration is 
considered, the shock velocities that were inferred from thermal line-widths 
may be grossly underestimated. Additionally, the fast shocks in Tycho, while 
not requiring a large cosmic ray population, do not exclude it either. If the 
outer blast-wave of Tycho is efficiently accelerating particles to cosmic ray 
energies, synchrotron radiation from high energy electrons may explain the 
featureless rim spectra while the thermal populations could be nearly equili¬ 
brated similarly to IE 0102—72. 

Progress in this developing held depends on accurate modeling of the emis¬ 
sion from pre-shock and post-shock gas, corroborating evidence from multiple 
wave-bands, and useful assessments of the cosmic ray production and its ef¬ 
fect on the shock. The spectral resolution of future X-ray missions such as 
ASTRO-E2, ConX and Xeus will allow thermal Doppler broadening to be 
measured routinely in the X-rays as they now are in the optical and ultravio¬ 
let. Observations of 500-2000 km s“^ shocks, where the transition from full to 
minimal equilibration should occur are especially important for diagnosing the 
relationship between shock velocity and the relative importance of different 
collisionless heating processes. 

I would like to thank Arend Sluis, Parviz Ghavamian and Kelly Korreck for 
fruitful discussions regarding this review and shock physics in general. This 
review also benehted from the comments of two conscientious referees. This 
work was supported by NASA grant NAG5-9281. 
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